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ABSTRACT 

In our previous paper, we evaluated the transit duration variation (TDV) effect for 
a co-aligned planet-moon system at an orbital inclination of i = 90°. Here, we will 
consider the effect for the more general case of i ^ 90° and an exomoon inclined 
from the planet-star plane by Euler rotation angles a, j3 and 7. We find that the TDV 
signal has two major components, one due to the velocity variation effect described 
in our first paper and one new component due to transit impact parameter variation. 
By evaluating the dominant terms, we find the two effects are additive for prograde 
exomoon orbits, and deductive for retrograde orbits. This asymmetry could allow 
for future determination of the orbital sense of motion. We re-evaluate the ratio of 
TDV and TTV effects, 77, in the more general case of an inclined planetary orbit with 
a circular orbiting moon and find that it is still possible to directly determine the 
moon's orbital separation from just the ratio of the two amplitudes, as first proposed 
in our previous paper. 

Key words: techniques: photometric — planets and satellites: general planetary 
systems — occultations — methods: analytical 



1 INTRODUCTION 

With exoplanet detection rates soaring, it is now becom- 
ing increasingly possible to characterise these alien worlds. 
Part of this characterisation will undoubtedly involve de- 
termining if exoplanets have their own moons, so called ex- 
omoons. The theoretical foundation s of exomoon detection 
using tr ansits were f ir st laid down bylSartoretti fe Schneider 
" 19991) . iDeesd l|2002t l. ISzabo et al l |2006l ) and ISimon et al. 
2007). In these papers, the emerging theme of using transit 
timing variation (TTV) as a detection tool was advocated 



In our previous paper. iKippind (|2009l ) (hereafter K09), we 
showed that an exomoon should induce not only a transit 
time variation (TTV) effect but also a transit duration vari- 
ation (TDV) effect on the host planet. The two effects were 
predicted to exhibit a 7r/2 phase difference which could be 
used as the hallmark signature of an exomoon. 

In our previous work, an underlying assumption was 
coplanarity within the system. We assumed that both the 
planet's orbital inclination angle, i, was 90° and that the 
moon's orbit was completely coplanar with the planet-star 
orbit. In this paper, we will extend the theoretical framework 
to include non-coplanarity. As a result of this consideration, 
we predict that the TDV effect due to an exomoon has infact 
two primary constituents: i) a velocity (V) component ii) a 
transit impact parameter (TIP) component. 

* E-mail: d.kipping@ucl.ac.uk 



The V-component is the same effect we described in our 
previous work, where the velocity of the planet is perturbed 
by the moon's presence. The TIP-component is a new effect 
which is due to the planet moving between higher and lower 
impact parameters as a result of the wobbling. Since transit 
duration is a strong function of impact parameter, then even 
slight changes can induce a TDV effect. 

This additional TDV component acts constructively 
with the V-component in the case of a prograde exomoon 
orbit and destructively for a retrograde orbit. With most 
large moons taking prograde orbits within the solar system, 
it would seem reasonable to expect constructive interference 
to be the typical case. As a result, the expected TDV signal 
from an exomoon is even more detectable. Furthermore, we 
predict this asymmetry could allow for a determination of 
the satellite's sense of orbital motion. 



2 THE TIP-COMPONENT OF THE TDV 
EFFECT 

In general, we posit that there exists two dominant com- 
ponents of the TDV effect. The V-component is caused by 
an exomoon perturbing the planet's velocity as it orbits the 
host star. This effect was discussed in depth in our previous 
paper and details can be found in K09. The second effect is 
the one we will concentrate on in this work and we label it 
as the transit impact parameter (TIP) component. 
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Figure 1. Cartoon of the side-on view of the star-planet-moon 
system. In this schematic, the star lies at the apex of the lines 
in the bottom left, the observer lies at y = +00 and the exomoon 
is not shown. The wobble of the planet is represented by the two 
grey spheres, (1) and (2), being the planet's maximal positions. 
The presence of a moon causes distance a perturbation in the 
distance q. 



Consider a side-on view of a planetary transit as shown 
in figu re 1. As demonstrated by ISeager fc Mallen-Ornelasl 
(|2003h . the transit duration is a strong function of the im- 
pact parameter of the transit, b, which is given by: 



6 =^ = t W 

The variable q is shown in figure 1 and denotes the 
distance between the observer's line-of-sight to the planet's 
centre, during the transit. It is this distance, q, which the 
transit duration is particularly sensitive to. If q — > 0, the 
planet transits across the star's equator, which is the star's 
widest point, thus giving a very long transit duration. If 
q — > R*, then the planet only grazes the star during the 
transit event and so we expect a very short duration. 

Now consider adding an exomoon. As discussed in our 
previous work, the exomoon is likely to be too small to ob- 
serve directly, but its gravitational effects on the host planet 
will be quite visible. Consider placing a moon around the 
planet such that the planet-moon orbital plane is the same 
as the star-planet orbital plane. Due to the moon's presence, 
the planet will exhibit reflex motion, or put simply a wobble, 
in this plane. From figure 1, it is clear that a component of 
this perturbation will be in the ^-direction. 

The motion along the 2-axis is of particular interest 
because we have already discussed how sensitive the transit 
duration is to the distance q. Any motion in this direction 
will cause q to get periodically larger and smaller. Ergo, the 
transit duration will vary. 

In the proposed case, we consider an exomoon with an 
orbit coaligned to a planet-star plane at an orbital inclina- 
tion angle i, where i ^ 90°. It is clear to see that there will, 
in general, always be a component of 2-axis wobbling motion 
for other moon inclinations. One of the special cases where 
this will not occur is when i — 90° and the moon takes a 
coplanar orbit, which is precisely the case we considered in 
our previous work. 



3 THE UPDATED MODEL 
3.1 Basic setup 

In our updated model, we consider a planetary orbit where 
i ^ 90°. We derive the total TDV effect by considering sev- 
eral stages of geometric manipulation of the planet's motion. 
We use the same reference axes in figure 1 where the observer 
is at y — +00, the planet-moon barycentre's true anomaly is 
fp, and the planet's true anomaly around the planet-moon 
bary centre is fw- In appendix A, we briefly consider the 
effects of inclined moon orbits, but for the mathematically 
simpler (and possibly more probable) case of a co-aligned 
moon orbit, it is shown that we may write the position of 
the planet as: 

z(fPi fw) = r P cos(/p + zu P ) + r w cos(f w + ww) 
z (fp,fw) = [r P sin(/p + zu P ) +rwsin(fw + ^w)} cos % 

(2) 

Where rp and rw are the planet to star and planet to 
planet-moon barycentre separations respectively. Note that 
we adopt the same notation as our previous papefl q is 
nominally given by rp cos i (see equation 1), but here we 
consider that q has been perturbed by the moon's presence 
to a new value given by q' . We denote the perturbation itself 
by Aq. It can be seen that q' = z(fp — > f m id,fw), where 
fmid is the true anomaly at the moment of mid-transit, given 
by fmid = tt/2 - tap- 

Without a moon, we would simply have z(fp) = 
rp cos i sin(/p + zup) and thus z(fp — > fmid) = q = rp cos i, 
i.e. no perturbation. Defining q = q + Aq, according to the 
directions in figure 1, we now have Aq — q — q: 

Aq — rw sin(/iv + ^w) cos i (3) 

Note that Aq > represents a shift in the +z-direction 
whereas Aq < is a shift in the — z-direction. As expected, 
for i — + 90°, we have Aq — > 0. For a highly inclined moon 
orbit, it is worth noting that the maximal value of Aq will 
be Aq — aw- 



3.2 Prograde versus retrograde orbital motion 

Consider the planet-moon barycentre moving in the +x- 
direction in figure 1. For a prograde orbit, the velocity of the 
planet around the planet-moon barycentre must be in the 
+a;-direction when it is at position (1). At position (1), the 
transit impact parameter has increased and thus the tran- 
sit duration has shortened. At the same time, the planet's 
wobble velocity is additive to the planet-moon barycentre 
velocity around the host star, and so the transit duration is 
further shortened. Thus for prograde orbits, it can be seen 
that the TIP- and V-components are additive. The opposite 
is true for retrograde orbits. 



1 Except for tilting the orbital plane into the +z direction rather 
than the negative 
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Table 1. Summary of key properties of the three known transit 
timing effects due to an exomoon. 



TTV 



TDV-V 



TDV-TIP 



Type of effect Positional 

Direction x 

Proportionality Msas 

Relative phase 



Velocity Positional 

x z 
1/2 



M s a s 
tt/2 



M s a s 
±tt/2 



3.3 Derived total TDV effect 

In appendix C, we evaluate the total TDV effect in the case 
of eg, a, (3, 7 = and ^ ep < 1. We are able to show that 
the r.m.s. amplitude of the TDV signal is given by: 



[ (R, 



2 • 

CLwCLp COS I 



2ixa 



(R, + R P ) 2 ~a 2 p cos 2 i 



T I P - com p o lie nt 



— (4) 
Ps v B± l ^2 



V-componcnt 



Where aw is the semi-major axis of the planet's orbit 
around the planet-moon barycentre, ap is the semi-major 
axis of the planet-moon barycentre's orbit around the host 
star, i is the orbital inclination angle of the planet-moon 
barycentre, vp± is the projected velocity of the planet-moon 
barycentre across the face of the star during transit, Ps is 
the orbital period of the satellite and f is transit duration 
of the planet in the absence of a moon. 

The positive sign refers to prograde moon orbits and the 
negative signs refers to retrograde orbits. It can be seen that 
the TDV effect has two dominant terms. The first term is 
the TIP-component is oc as Ms, which is the same as TTV's 
proportionality (where Ms is the mass of the satellite). The 

— 1/2 

second term is the V-component and is oc Msa s , as found 
in our previous paper. A summary of the properties of the 
three known transit timing effects due to an exomoon can 
be seen in table 1. 




0.4 0.6 0.8 

Transit impact parameter, b 

Figure 2. For a hypothetical planet-moon system, we plot the 
variation of the TDV effect with transit impact parameter. Whilst 
partially transiting planets seem to offer enhanced signals, they 
are also highly unlikely to be detected and thus the case of b = 
is more likely to be the optimal condition. 



In figure 2, we plot the variation of the TDV effect for 
a prograde 1M$ exomoon around a hypothetical planet. We 
use the same hypothetical planet as in our previous, K09, i.e. 
an identical system to GJ436b except the orbital period is 
35.7 days and the eccentricity is zero. As seen in figure 2, the 
TDV effect slowly drops off for increasing b until we reach 
the partial-transit regime where the TDV effect becomes 
extremely large due to the TIP-component dominating. 

However, the probability of detecting a transit which is 
only partially transiting is very small due to two reasons: 1) 
geometrically the inclination range is very small 2) V-shaped 
lightcurves are usually rejected as a planetary candidate and 
labelled as a grazing eclipsing binary. For these reasons, we 
consider the optimal planets for detecting an exomoon to 
have the longest transit duration possible, i.e. b ~ 0. 



4 IMPLICATIONS 

4.1 TDV's inclination dependence 

By considering orbital inclination and the transit impact 
parameter, we have demonstrated that an additional TDV 
effect exists, which we have labelled as the TIP-component. 
It may be tempting to assume that inclined orbits therefore 
improve the TDV effect due to an exomoon (assuming a pro- 
grade orbit) and perhaps the best exoplanets candidates for 
detection purposes would be near-grazing transits. However, 
we counter that naive supposition by pointing out that the 
transit duration itself decreases with higher impact parame- 
ters and it can be seen in equation (4) that both the V- and 
TIP-components have dependancies on f(b). 

We may rightly ask whether there is a certain value of 
impact parameter which enhances the TDV signal optimally. 
By differentiating equation (4) with respect to b and making 
some approximations, we find that the optimal value of b 
occurs for b ~ (I + k), where k is the ratio-of-radii. In other 
words, the TDV signal is enhanced for partial transits where 
the impact parameter is so high the lightcurve takes a V- 
shape. 



4.2 Determining Ps and Ms 

What are the consequences for determing both as and Ms 
separately, as proposed in our last paper? Certainly taking 
the simple ratio of the TDV and TTV effect will not provide 
precisely the s ame equation as we predicted before. From 
iKippingj l|2009l ). the TTV r.m.s. amplitude, for a circular 
moon orbit, is given by: 

. I asMsJap 1 ._. 
OTTV ~ (5) 

\/2 M P Rvy/G{M* + Mprv) T 

Where T accounts for the planet's orbital eccentricity. 
Taking the ratio of TDV to TTV, we have: 

_ Stdv 
Sttv 

_ y/Ga P (M* + Mp RV )Tf cos 2 i 2ttt 

11 ~ {R* + R P ) 2 - a 2 p cos 2 i ^Ps^ [ ' 

'constant' 'info' 

The TIP-component makes r\ change from simply being 
inversely proportional to Ps to being inversely proportional 
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Table 2. Predicted TTV and TDV (both V- & TIP- components) rms amplitudes due to a LMgj exomoon at 1/3 the Hill radius, for a 
selection of the best candidate transiting planets. System parameters are taken from various references, which are shown. 



Planet 


Sttv I s 


V-part of &tdv/s 


TIP-part of 


HAT-P-llb 


19.19 


22.54 


0.40 


GJ436b 


14.12 


13.68 


1.30 


CoRoT-Exo-4b 


7.58 


9.15 


0.00 


OGLE-TR-lllb 


4.63 


7.32 


0.11 


HAT-P-lb 


4.58 


6.82 


0.47 


HD149026b 


3.61 


9.76 


0.00 


Lupus-TR-3b 


3.28 


5.19 


0.07 


WASP- 7b 


3.26 


5.88 


0.00 


HD17156b 


3.07 


1.06 


0.43 


TrES-lb 


3.04 


5.95 


0.05 


HD209458b 


2.97 


5.95 


0.07 


XO-5b 


2.65 


4.69 


0.17 


HAT-P-4b 


2.54 


8.34 


0.00 


HD189733b 


1.52 


2.96 


0.16 


XO-3b 


0.41 


0.87 


0.07 



rl* 



Reference 



Bakos et al. 



Alonso et al. (2008) fe Torres et al. (2007) 
Aigrain et al. (20081 

Diaz et al. (2008) 
Johnson et al. (2008) 

Winn et al. (20071 
Weldrake et al. (2008) 
Hellier et al. (2009) 
Barbieri et al. (20071 
Winn et al. (2007a) 

Kippine (_2008) 
Burke et al. (20081 
Kovacs et al. (20071 
Winn et al. (2007b 1 & Beaulieu et al. (20081 
" Winn et al. (20081 



to Ps plus a constant. We highlight the following key points 
about this equation: 



■ The 'constant' is not a function of the moon's properties 
whatsoever, i.e. it depends on the planetary properties only. 

■ In the case of i = 90°, the 'constant' — > 0, retrieving 
the original equation for r), as derived in K09. 

■ The 'constant' may be calculated independently of the 
moon's properties in a reliable way and thus Ps may still be 
calculated, providing we assume es = and a, j3, 7 = 0. 



Equation (6) tells us that it is still possible to evaluate 
P s by simply taking the ratio of the TDV and TTV effects, 
provided we make certain assumptions about the moon's 
orbit. The 'constant' term quoted above is not a function 
of the moon's properties, but in reality it is actually a very 
weak function since we have made the assumption ap> aw , 
which is certainly a very valid assumption to make. 

Although a large TIP-component can enhance the de- 
tectable signature of an exomoon, if the TIP-component is 
greatly larger than the V-component, then our ability to ac- 
curately estimate Ps and hence Ms will diminish. Ideally, 
the V-component should remain the dominant term for ac- 
curate determination of these parameters. 

For a system of interest like GJ436b, the constant is 
~ 0.1 suggesting that TTV is an order of magnitude stronger 
than the TIP-component of the TDV effect. Although this 
is small, the key point is that it is additive and can be larger 
than the V-component in certain special cases, notably graz- 
ing transits. 

We conclude that our previous statement that the ra- 
tio of TDV and TTV can be used to determine the moon's 
orbital distance and mass still remains true, provided we as- 
sume a circular co-aligned exomoon orbit. For systems with 
exomoon eccentricity and inclinations, there will be insuf- 
ficient information through timing alone to solve for all of 
these parameters. 



4.3 Determining the sense of an exomoon's 
orbital motion 

Observations of transit timing may also permit the determi- 
nation of whether an exomoon is in a prograde or retrograde 
orbit, given sufficient signal-to-noise. We will illustrate this 
possibility by referring to a hypothetical example of an ex- 
omoon detection. We use the same example as in our last 
paper, K09, where we considered shifting GJ436b to an cir- 
cular orbit of 35.7 days period and add a 1M© exomoon 
on a 2.5 day orbit. We use the same impac t parameter as 
measured for GJ436b bv lTorres etafl (|2007h . i.e. b = 0.848. 

Consider we measure the TTV r.m.s. amplitude of this 
planet to be 137.4±0.4 seconds and the total r.m.s. TDV am- 
plitude would be 39.7±0.8 seconds giving r) = 0.289± 0.020. 
Note that we have assumed the duration is measured to half 
the precision of the mid-transit time and assumed timing 
errors in-line with the capabilities of forthcoming missions. 
Based on the known inclination and planetary properties, we 
are able to evaluate the r\ 'constant' term to be 0.02739 to 
a negligible error (since this can be refined by compositing 
multiple transits). Therefore the 'info' component of r\ is ei- 
ther 0.261±0.020, if prograde, or 0.316±0.020, if retrograde. 
This corresponds to an exomoon period of either 2.52 ±0.14 
days or 2.08 ±0.14 days for prograde and retrograde respec- 
tively, differing by 3-sigmas. With the TTV measured to a 
signal-to-noise of over 350, it would not be difficult to use 
multiple TTV measurements to identify which of these peri- 
ods is permitted by the frequency of the data points, which 
in this case is the prograde orbit. 

It is important to remember that this calculation was 
done for a planet-moon orbital plane which is coaligned to 
the star-planet orbital plane. Slight moon inclinations of 
< 5° would not change the result significantly but very 
large inclinations would severely disrupt this calculation's 
accuracy. We propound that exomoons of low inclination 
angles would be identifiable by a planet-moon eclipse which 
should be observable in t he lightcurve, as predicted by 
ICabrera fc Schneider! (|2007l ). 

We therefore propose that it should be possible for fu- 
ture observations to not only detect an exomoon and deter- 
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mine its mass, but also provide a confident deduction of the 
sense of orbital motion. Although this determination will 
likely require photometry at the limit of planned missions, 
it seems likely that once an exomoon is detected a more in- 
depth investigation would be able to answer the question of 
orbital sense of motion conclusively. 



5 DISCUSSION AND CONCLUSIONS 

We have shown that an exomoon around a transiting exo- 
planet should induce a transit duration variation effect with 
two dominant components. One of these components is due 
to the moon altering the velocity of the host planet, which 
we label as the V-component. The second constituent is due 
to the impact parameter of the transiting planet varying 
as a result of the moon's presence, which we label as the 
TIP-component . 

In table 2, we have evaluated the V- and TIP- com- 
ponents for a list of targets as in our previous paper. The 
table suggests that the TIP-component is often an order 
of magnitude less than the V-component, but can therefore 
we several seconds for some targets. We do not anticipate 
this additional component to be a hurdle in determining the 
moon's mass and orbital distance since the dominant effect 
on the rj parameter is to introduce an additive constant, 
which is independent of the moon's properties. 

The TDV effect can be markedly increased for prograde 
moons which improves their detectability. For grazing tran- 
sits with highly inclined moons, the planet could even go 
through epochs where it no longer transits at all, although 
we do not expect this to be a typical situation. Further- 
more, we predict that including the TIP-component may 
allow for the determination of the orbital sense of motion 
of an exomoon around an exoplanet. We do however stress 
that such a determination would require very precise, but 
feasible, photometry. We believe this paper further demon- 
strates the feasibility of detecting an exomoon, but outlines 
the great care and understanding required to complete the 
analysis. 



ACKNOWLEDGMENTS 

DMK is supported by STFC, University College London and 
HOLMES ANR-06-BLAN-0416. The author would like to 
thank Jean Schneider for technical discussions on the sub- 
ject of transit timing effects. Special thanks to the reviewer 
of this paper whose comments have greatly improved the 
manuscript . 



Cabrera, J. & Schneider, A&A, 464, 1133 
Deeg, H. J., 2002, ESA SP-514, 237 

Diaz, R. F., Rojo, P., Melita, M, Hoyer, S., Minniti, D., 

Mauas, P. J. D., Ruiz, M. T. 2008, ApJ, 682, L49 
Domingos, R. O, Winter, O. C, Yokoyama, T. 2006, MN- 

RAS, 373, 1227 
Hellier, C. et al. 2009, ApJ, 690, L89 
Johnson, J. A. et al. 2008, ApJ, 686, 649 
Kipping, D. M., 2008, MNRAS, 389, 1383 
Kipping, D. M., 2009, MNRAS, 392, 181 (K09) 
Kovacs, G. et al. 2007, ApJ, 670, L41 
Sartoretti, P. & Schneider, J., 1999, A&AS, 14, 550 
Seager, S. & Mallen-Ornelas, G, 2003, ApJ, 585, 1083 
Simon, A., Szatmary, K., Szabo, Gy. M., 2007, A&A, 470, 
727S 

Szabo, Gy. M., Szatmary, K., Diveki, Zs. & Simon, A., 2006, 

A&A, 450, 395 
Torres, G, et al., 2007, ApJ, 666, L121 
Weldrake, D. T. F. et al. 2008, ApJ, 675, L37 
Winn, J. N., Henry, G. W., Torres, G & Holman, M. J. 

2007, ApJ, 657, 1098 
Winn, J. N., Holman, M. J. & Roussanova, A. 2007, ApJ, 

657, 1098 

Winn, J. N., Holman, M. J., Henry, G. W. et al. 2007, AJ, 
134, 1707 

Winn, J. N. et al. 2008, ApJ, 683, 1076 



APPENDIX A: MATHEMATICAL 
TREATMENT 

In appendix A, we will derive the two components of the 
TDV effect separately. The two effects will later be combined 
in appendix B. The total TDV effect will then be evaluated 
in the case of a circular, coaligned moon in appendix C. The 
error regarding one of our key assumptions will be calculated 
in appendix D. 

Al Planetary wobble motion 

From figure 1, it is clear that the change we care about, 
Aq is in the +z- direct ion. In order to create a generally 
oriented orbital plane, we do so in several steps. In our first 
step, we consider the planet-moon centre of mass frame and 
employ the same geometric model as that of our previous 
paper where we consider an elliptical orbit with a centre at 
the origin of an x'-y'-z' co-ordinate system, the S' frame. 
We may write the position of the planet, as a function of its 
true anomaly within this frame (fw), as: 
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where 



x' = awe-w + rw cos fw 
y = r w sin fw 
z'=0 



rw 



a w (l — e w ) 



(Al) 



(A2) 



1 + ew cos fw 

In the same manner described by iKippin j (2008), 
we then transform these positional coordinates with i) a 
counter-clockwise rotation about z'-axis by the position of 
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Figure Al. Euler angles of our rotational scheme. 



pericentre angle, -cow (the S[ frame) ii) a translation of the 
planet-moon barycentre to the origin (the 5*2 frame). These 
two transformations give us: 

x 2 = rw cos(/w + zuw) 
y' 2 = r w sin(/W + tew) 
4=0 (A3) 

It is worth noting that in the case of ew — 0, which we 
would anticipate to be the typical scenario, we have x' 2 = 
aw cos fw and y' 2 = aw sin fw- 

In this system, we allow the orbital arrangement of the 
moon to be rotated by the three Euler angles, a, j3 and 7 
with respect to the star-planet plane. We choose to employ 
the conventional z-x-Z convention for the Euler system. The 
Euler rotation angles are ideal because any general rotation 
in three dimensions can always be written in terms of these 
three angles. 

We denote the x' 2 ,y' 2 ,z 2 position of the planet away 
from the planet-moon barycentre origin by the vector X 2 , 
then after the Euler rotations we have a new positional vec- 
tor given by: 

X 3 = R Eulcr (a,/3, 7 )X 2 (A4) 

where R.Euicr is the Euler rotation matrix. Since z' 2 = 0, 
we may write: 

X3 — [c a C 7 — S a Cf3Sj] ■ X'2 + [— C a S 1 — SaCpCj] ■ y' 2 
y'z = [SqC 7 + C a Cf3S~,] ■ X 2 + [—S a Sj + C a C0C~,] ■ y' 2 

z' z = sps~ f ■ x 2 + Sgc~, ■ y 2 (A5) 



We have therefore found the position of the planet away 
from the planet- moon barycentre as a function of fw- 

A2 Planet-moon barycentric motion 

Consider the orbit of a point mass around the host star. 
We denote this point mass by the coordinates x, y, z (the S 
frame). Placing the origin at the centre of the ellipse, we can 
write: 

x — apep + rp cos fp 
y = r P sin f P 

5 = (A6) 

We then rotate for position of pericentre and translate 
so that the star is at the origin. 

x 2 — rp cos(/p + vjp) 
y~2 = r P sin(/p + zu P ) 
£2 = (A7) 

For the moment we will not include the orbital inclina- 
tion angle, but return to it later. 

A3 Overall motion of planet 

Earlier on we defined the position of the planet, in relation to 
the planet-moon barycentre as x' 3 ,y' 3 ,z' 3 . So the true position 
of the planet is the vector for the barycentre's position, X2 
added to the vector for the relative position of the planet, 
with respect to the barycentre, X3. This gives us Xa = 
X2 + X 3 . 

XA = X2 + x' 3 

yA = 2/2 + y' 3 

z A = zi + z 3 (A8) 

Finally, we now rotate about the iA-axis in a counter- 
clockwise sense by an angle (it/ 2 — i), in order to be consis- 
tent with our defined system shown in figure 1. Note that in 
our previous paper we performed this rotation in the clock- 
wise direction and we choose to do the opposite direction 
here for mathematical simplicity. Of course, the choice of 
direction does not affect the final result. 

X = X A 

y — yA sin i — za cos i 

z — yA cos i + za sin i ( A9) 

A4 Application to a co-aligned moon orbit 

Although the Euler angles can play a significant role in en- 
hancing the TDV effects, we predict that most exomoons 
will not exhibit significant inclination deviations from the 
planet-star plane and so we decide to set all of these angles 
to zero. This is also valid because it would be impossible 
to determine all of the Euler angles from just two measure- 
ments in any case. Thus we have: 
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X — x 2 + x 2 

y = [m + y's] sini 

z = [y~2 + y' 2 ] cos i 



(A10) 



For an observer at infinite +y, the {/-component is never 
seen. So we have: 



x (fp, fw) = r P cos(/p + ujp) + r w cos(f w + vcw) 

z (fp,fw) = [rp sin(/p + vj p ) + r w sm{f w + ncrw)] cosi 

(All) 

From figure 1, q' = z(f P -» f mid , f w ) and Aq = q' - q 
where q is the offset in the z-direction when no moon is 
present. Removing the perturbation terms, we can derive 
the moonless quantity to be q — rp sin(f m id + top) cos i = 
rp cos i. The two equations therefore give the perturbation 
magnitude to be: 

Aq — rw sin(/w + vjw) cos i (A12) 



A5 TIP-component of the TDV effect 

Here, we will derive approximate equations for the TIP- 
component of the TDV effect. For simplicity, we may write 
down the case of ew = 0. We also assume a = f3 — 7 = 0, 
but i 7^ 90°. Ergo, this is not the same case considered in our 
previous research where i = 90°. Although we could write 
the down the expressions for the most general case, to do 
seems pointless given that there will not be enough informa- 
tion to solve for all of these parameters in any case. For the 
stated assumptions: 



Aq(fw) = aw cos i sin fw 
vw±(fw) = ±|«wJ-| sin/iv 



(A13) 



where the ± symbol refers to prograde / 'retrograde exo- 
moon orbits respectively 

Although in reality the magnitude of the vector con- 
necting the planet to the star's centre is no longer ap, the 
maximum by which it can be modified will be given by aw- 
Since ap 3> aw, we will assume that this changing magni- 
tude vector does not significantly affect the TDV amplitude 
relative to the change in q. 

As a result of changing q, the observed transit impact 
parameter appears to change. Without a moon present, the 
impact parameter should be simply given by: 



q 



(A14) 



But now the altered effective impact parameter will be 
given by: 



b'(fw) = ^- = 



R* 



(A15) 



We now consider that the TDV signal due to the TIP 
effect is given by the observed duration (r) minus the ex- 
pected duration: 



TDV(jW) = rift) - r(6) = r{b') 



(A16) 



A6 Modification to the V-component 

iKippin j ((2009) proposed that the velocity of a planet during 
a transit is modified by the presence of an exomoon. We 
label this TDV effect as the V-component. This effect is 
in addition to the previously detailed TIP-component. We 
will now derive the modification to the V-component in the 
presence of orbital inclinations. In our previous paper, we 
defined: 



TDV(fw) = r(f w ) 
lim TDV(fw) = 



VB± 



(A17) 



*t/2 ' \V B ± + VwJ-(fw) 

In our orbital setup, we assume the planet is moving in 
the positive x-direction. Therefore, we need to know what 
component of vw remains in this direction given a rota- 
tion through the three Euler angles. Taking the original 
vector vw = {r>w±,0,0} and rotating we find that the x- 
component is modified to: 

r>w± — * [cos (a) cos(7) — sin(a) cos(/3) sin(7)]vv7i 
vw± — » <p(a, (3, 7) • v w ± (A18) 
This modifies our velocity TDV component to: 



lim TDVf/W) 

•;— >7r/2 



iifli 



+ tp(a,P,i) ■ vw±(t) 



VB- 



(A19) 

But note that this does not include the TIP-component. 
Now that we have written down the TDV signal for both the 
V- and TIP-components independently, the next step is to 
combine the two. 



APPENDIX B: THE TOTAL TDV EFFECT 

We now consider the total effect. The V-component is just 
a factor which modifies the duration, so this can be applied 
after the TIP-component. By employing this ordering, we 
are able to write the transit duration, in the the general 
case, as: 



VB± 



vb± + p(a,(3,j) ■ v w ±(fw) 



) ■ r(6') (Bl) 



It would useful at this point to have r(b') written as 
some factor multiplied by the mean transit duration f , where 
it is understood that b'(fw)- Let: 



e(b') 



(B2) 



We note that r gets larger as b approaches zero. This 
occurs for Aq being positive. Thus, e > 1 for Aq > and 
vice versa. Our TDV may now be written as: 



TDV(/ W 



s(b') ■ v B ± 



v B ± + <p(a,/?,7) ■ vw±(fw) 



(B3) 
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Appreciating that vw± -C Vp±, 
approximated to: 



this expression may be 



TDV(/ W 



1) 



£ ■ Vw± 

"SI 



(B4) 



We can now see that the TDV signal has two clear com- 
ponents. In the absence of any TIP-component, e — > 1 and 
hence we recover the original TDV effect predicted in our 
previous paper. Since e is a value close to unity, we choose 
to write it as: 



1 + Q 



(B5) 



Where it is understood that g is small compared to 
unity. This now gives us: 



TDVCM- [g- 
TDV(.M * [g(fi 



V\V± 



VB± V B ± 

n vw±(fw 



VB± 



(B6) 



Consider the case of Aq > 0, as mentioned earlier this 
means g > and so the first term is positive. Aq is posi- 
tive when vw{fw) goes negative for prograde orbits. So the 
second term must be negative too, and thus we have a dou- 
ble negative which equals a positive. Thus we confirm that 
when for prograde orbits the TDV effect is additive. 

It is also clear that the total effect is dominated by two 
terms which are linearly additive. These means the TIP- 
component can significantly increase TDV signals due to 
prograde exomoons. On the other hand, it can significantly 
dampen any effect for retrograde moons. Infact, if the TIP- 
component is very large for a retrograde, it may change the 
sign of the effect completely meaning that instead of TDV 
lagging behind TTV by it/2, the opposite is true. 



e{b'f = 



1 + ^) -Ok) -2(^)-0(^ 



R- 



R* 



,. \ 2 / \ 2 
Rp \ I _q_ 

R, 



(C3) 



write: 



Assuming O^j^j is small and ap S> aw, we ma y 



e(b') 



2awap cos 2 i sin fy 
(R*+R P ) 2 - a 2 , cos 



(C4) 



Unfortunately, evaluating the integral of TDV 2 is non- 
trivial, even with these approximations. Without a direct 
computation of this integral, we cannot evaluate the r.m.s. 
amplitude directly. In order to estimate the r.m.s. ampli- 
tude, we assume the signal takes a sinusoidal form and thus 
the r.m.s. amplitude will be given by the normal amplitude 
divided by y2~. This approximation is particularly valid be- 
cause we have assumed es = 0. The maximum value of e and 
vw± is for fw = — 7r/2. We also take advantage of the fact 
we know e is close to unity and therefore we approximate 
the square root to find g ma x- 



1 + 



2 • 

aw a p cos i 



Qrnax 



(R, + R P ) 2 -a 2 p cos 2 



aw a p cos i 



— 1 ~\~ Qmax 



(R*+Rp) 2 - a 2 , cos 2 i 
This gives us a max signal of: 



(C5) 



TDV 



l(R* 



aw ap cos i 



2ixa\ 



[R, + Rp) 2 
Or an r.m.s. signal of: 



a p cos" i 



PS Dfll 



— 1 r (C6) 

3± J 



APPENDIX C: 
EFFECT 



EVALUATION OF THE TDV 



CI Case of e s ,a,/3,j = & < e P < 1 

In order to ascertain the proportionality and magnitude 
of the total TDV effect, let us first assume the transit 
duration is given by the circu lar equations described by 
ISeager fc Mallen-Ornelaii (|2003T ). Using equation (16) from 
this paper and replacing ap cos i with q: 



P P R, 

T ~ 

nap 
Which gives us: 



1 + 



Rp\'< 
R* ) 



(CI) 



+ £ 

'1 + 



Rf 



2 / \ 2 
1 

R„. 



(C2) 



Even if we use the equations for the approximate transit 
duration due to an eccentric orbit as presented by ? by their 
equation (1), the same value of s 2 is derived. Thus the fol- 
lowing derivation holds true for planets on eccentric orbits. 
Expanding out the brackets we have: 



5td\ 



awap cos i 



2-Ka\ 



.(R,+R P ) 2 



Ps v B ± J ^2 



(C7) 



APPENDIX D: ERROR ON CONSTANT 
VELOCITY ASSUMPTION 

Throughout we have assumed that vw± does not change 
significantly over the course of the transit. This essentially 
the same as assuming that Ps 2> f. For a normal circular 
orbit, the planetary transit would occur for fp changing by 
a quantity: 



Afp — arcsin |^ 



(l + fc) 2 + (a P cosi/7?*) : 



(Dl) 



If the moon and planet are on a circular orbit, and Ps = 
HPp, then the change in fw may be written as: 



a f 1 . r(l + fc) 2 + (a P cosi/i?*) 2 l 

Ajw = — ■ arcsin ^5 ■ — (D2) 

a L sin i J 

The velocity of the planet wobble is given by: 



vw± = \v w ±\ sin fw 



(D3) 



So the average velocity over the course of Afw is given 

by: 

. \VW±\ ffw,mi d +0.^fw 

< v w ± > (fw,mid) = . . / sm/w Afw 

^J W J fw,mid-0-5Afw 

jsin(A/ w /2)- 



r sm(A/iy/2) i 
|^x|sm( /tv ,^)[ A/W2 J 



(D4) 

In contrast, previously we assume a fixed constant ve- 
locity of vw± = \vw± | sin fw,mid- So it seems we very 
slightly overestimate the average velocity during the transit, 
as expected. The dominant error term can be written out by 
expanding the sin function: 



Vw±(fw) — \vw± \ sin(f w ,mid) — \vw±\ sin(/ lv , mic j 



Afw 



(D5) 

Thus the fractional error in Vw± is given by (A/^/48). 
The error should be largest for close-in orbits and even for a 
3-day hot- Jupiter with typical parameters we would expect 
this error to be less than 1 part in 10 3 . 
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